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ABSTRACT 

With recent and archival Rossi X-Ray Timing Explorer (RXTE) X-ray measurements of the heavily 
obscured X-ray pulsar EXO 1722-363 (IGR J17252-3616), we carried out a pulse timing analysis to 
determine the orbital solution for the first time. Using a single pulse period derivative, we connected 
datasets separated by over 2.4 years without pulse number ambiguity (> 180,000 pulses), implying 
continuous spin-up of the pulsar from 2003 August (and possibly earlier) until 2006 February (and 
possibly later). The orbital solution also shows that a torque reversal occurred sometime between 
1998 November and 2003 August. The binary system is characterized by a^sinz = 101 ± 3 lt-s and 
-Porb = 9.7403 ±0.0004 days (90% confidence), with the precision of the orbital period being obtained 
by connecting datasets separated by more than 7 years (272 orbital cycles). The orbit is consistent 
with circular, and e < 0.19 at the 90% confidence level. The mass function is 11.7 ± 1.2 Mq and 
confirms that this source is a High Mass X-ray Binary (HMXB) system. The orbital period, along 
with the previously known ~414 s pulse period, places this system in the part of the Corbet diagram 
populated by supergiant wind accretors. Using previous eclipse time measurements by Corbet et al. 
and our orbital solution, combined with the assumption that the primary underfills its Roche lobe, we 
find i > 61° at the 99% confidence level, the radius of the primary is between 21 Rq and 37 Rq, and 
its mass is less than about 22 Mq. The acceptable range of radius and mass shows that the primary is 
probably a supergiant of spectral type B0I-B5I. Photometric measurements of its likely counterpart 
are consistent with the spectral type and luminosity if the distance to the system is between 5.3 kpc 
and 8.7 kpc. Spectral analysis of the pulsar as a function of orbital phase reveals an evolution of the 
hydrogen column density suggestive of dense filaments of gas in the downstream wake of the pulsar, 
with higher levels of absorption seen at orbital phases 0.5-1.0, as well as a variable Fe Ka line. 
Subject headings: X-rays: binaries — pulsars: individual (EXO 1722-363=IGR J17252-3616) 



1. INTRODUCTION 

EXO 1722-363 (IGR J17252-3616) was first discovered 
by EXOSA T scan observations in 1984 along with 7 other 
sources in the Galactic ridge and 14 sources outside the 
Galactic ridge (Warwick et al. 1988). A pointed Ginga 
observation in 1987 October for a total of 80 ks discov- 
ered pulsations with a 413.9 ± 0.2 s period (Tawara et al. 
1989). The source exhibited significant variability, with 
the 6-21 keV flux decreasing from about 2 mCrab to 
0.2-0.3 mCrab over 8 hrs. Using an absorbed power-law 
model, they found TVh = 1.2615^ x 10 24 cm~ 2 , suggest- 
ing a dense envelope of circumstellar matter surround- 
ing the pulsar. Additional Ginga observations were per- 
formed in 1988 March-April by Takeuchi et al. (1990). 
Using pulse timing analysis, lower limits on the orbital 
period and mass of the companion were set to 9 days and 
15 Mq, assuming a canonical neutron star mass of 1.4 
Mq. Takeuchi et al. (1990) also measured the spectrum 
but included a high energy cut-off, and found a lower 
column density N H = 1.26±g;|| x 10 23 cm~ 2 . 

Using RXTE Proportional Counter Array (PC A) scans 
of the Galactic center region from 1999 February to 2003 
October, Corbet el al. (2005) measured a 9.741 ± 0.004 
day orbital period through periodic changes in the source 
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flux. By folding the scan observations on this orbital pe- 
riod, eclipses were detected lasting 1.7 ± 0.1 days (eclipse 
half-angle 6 e = 31?8 ± 1?8). Corbet et al. (2005) also 
measured the spectrum using pointed PCA observations 
from 1998 October-November, and found that the spec- 
trum varies but not dramatically, although the spectrum 
was sampled over limited ranges of orbital phase. 

EXO 1722-363 has also been observed by INTEGRAL 
during surveys of the Galactic plane. Zurita Heras et 
al. (2006) found that the source is persistent in the 20- 
60 keV band but is not detected above 60 keV. On a 
time scale of ~3 days, the 20-60 keV source flux varied 
by a factor of about 4. Flaring activity was also ob- 
served, with the source flux increasing by a factor of 5 
and returning to its base level in about a day. Using 
INTEGRAL, the pulse period was measured to be 413.7 
± 0.3 s in 2003 August. Thus, aside from slight (~1 s) 
changes, the neutron star pulse period has been stable 
for over a decade. Phase-resolved spectroscopy did not 
show any variation of the continuum emission with pulse 
phase. Using a follow-up XMM-Newton observation from 
2004 March, a 414.8 ± 0.5 s pulse period was measured. 
More importantly, however, XMM-Newton measured the 
source position to 4" accuracy, and an IR counterpart 
2MASS J17251139-3616575 is located 1" away from the 
best position. The possible IR counterpart is not de- 
tected in the 2MASS survey in the J-band but appears 
in the -ff-band and JT s -band with 11.8 mag and 10.7 mag 
respectively. 

In this work, we use recent and archival RXTE PCA 
observations to determine the orbital parameters accu- 
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rately for EXO 1722-363 using pulse arrival time analy- 
sis. Since the orbital period for this system was previ- 
ously known, the new observations were timed to provide 
even coverage in orbital phase, allowing us to obtain firm 
constraints on the orbital parameters, and especially the 
eccentricity. The observations are described in § 2, as 
well as the pulse arrival time measurements and error 
analysis. The results of the pulse timing analysis are pre- 
sented in § 3.1. We confirm the 9.74 day orbital period 
and increase its precision by an order of magnitude. The 
orbital solution is then used to observe orbital modula- 
tion of the flux, pulse fraction, and pulse profile in § 3.2, 
to constrain the mass and radius of the mass-donating 
primary in § 3.3, and to study the orbital modulation of 
the source spectrum in § 3.4. We discuss the results of 
the spectral analysis in § 4. Finally, in § 5 we summarize 
our main results. 

2. OBSERVATIONS 

To find the orbital solution to EXO 1722-363, we ob- 
tained twenty new RXTE observations (ObsID 91080- 
02) at carefully chosen intervals during 2006 January to 
March (hereafter referred to as epoch 3 4 ), with each ob- 
servation roughly 3-5 ks in duration. In addition, we 
used two archival observations: One archival dataset 
(ObsID 30142-01) contains five observations from 1998 
October to November at ~4 day intervals (epoch 1), 
while the other (ObsID 80424-01) contains two observa- 
tions from 2003 August of 1 and 2 ks, respectively (epoch 
2). A summary of these observations is shown in Table 
1. 

Source and background light curves were created with 
data from the PCA (Jahoda et al. 2006) using stan- 
dard FTOOLS. Due to the faintness of the source and 
short observations, HEXTE data were not included in 
the analysis. The PCA instrument consists of five iden- 
tical multianode proportional counter units (PCUs), op- 
erating in the 2-60 keV range, with an effective area of 
approximately 6500 cm 2 and a 1° FWHM field of view. 
Photon energies for pulse timing analysis were restricted 
to 3-24 keV, and the light curves were binned by 1 s. 
Data accumulated in the Good Xenon data mode was 
used when available. The Faint 5 background model was 
used for the all epochs. The arrival times of individual 
events were reduced to the solar system barycenter using 
the Jet Propulsion Laboratory DE-200 ephemeris and the 
FTOOL faxbary. Figure 1 shows representative epoch 
3 light curves at three different points in the orbit. Peri- 
odic decreases in source flux, called dips, are observed in 
sufficiently inclined systems that have an accretion disk. 
These periodic decreases are thought to be caused by 
a bulge of material where the accreted matter impacts 
the disk, obscuring the X-ray source. Such dips typically 
last for several hours. By examining the low-energy light 
curve in detail (where photoelectric absorption has the 
largest effect), in addition to the hardness ratio, we do 
not find any evidence for dips in EXO 1722-363. 

2.1. Arrival Times 

Pulse arrival times were measured in a multi-step it- 
erative process, similar to the analysis in Thompson et 

4 Not to be confused with PCA calibration epochs. 

5 See http://heasarc.nasa.gov/docs/xte/pca_news.html 
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Fig. 1. — Representative light curves at three points in the 
pulsar's orbit for 3—24 keV photons using 8 s bins. The orbital 
phases are inferred from the orbital solution in § 3.1. 



al. (2006). First, each individual observation was di- 
vided into ~1.5 ks "light segments" (~3-4 pulse cycles), 
and each segment was folded modulo 414.0 s to improve 
statistics. Although the true pulse period may differ 
slightly from this value, we expect smearing of the pul- 
sation to be minimal, leading to systematic errors in the 
measured pulse arrival times of < 1 s. EXO 1722-363 is 
known to be an eclipsing system (Corbet et al. 2005). 
Zurita Heras et al. (2006) found that the INTEGRAL 
ISGRI flux drops to zero during the eclipse. The PCA 
flux, on the other hand, does not drop to zero during 
eclipse, probably because of Galactic ridge emission or 
possible weak contaminating sources within the 1° field 
of view. To limit the potential for including spurious 
pulse minima when the pulsar is obscured by the high- 
mass companion, each observation was fit with a constant 
flux and \ 2 was calculated. Based on the null hypoth- 
esis that the light segments are consistent with a con- 
stant flux, pulse minima were only measured if the null 
hypothesis was rejected at the 1% level of significance. 
While this may exclude data when the pulsar is out of 
eclipse, we prefer to be conservative and potentially ex- 
clude usable segments rather than include segments in 
eclipse. With the orbital solution obtained in § 3.1, we 
find that all excluded light curves are within |<IVb| < 0.15, 
where mid-eclipse is defined to be orbital phase 0.0. The 
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TABLE 1 

RXTE Observation Log of EXO 1722-363 



ObsID a 


Epoch 


Date 


Exp. Time 


Count Rate b 


Pulse Min. c 


Orb. Phase d 


Spec. Num.° 






(U.T.) 


(ks) 


(Counts/s/PCU) 






30142-01-02-01 


1 


1998 Oct 23.08-23.25 


14.80 


24.77 


10 


0.67 




30142-01-02-00 


1 


1998 Oct 27.14-27.31 


15.14 


6.30 


1 


0.09 




30142-01-03-00 


1 


1998 Oct 31. 93-01. 10 f 


14.22 


13.13 


7 


0.58 




30142-01-04-00 


1 


1998 Nov 05.00-05.17 


14.16 


4.18 





0.00 




30142-01-05-00 


1 


1998 Nov 09.07-09.25 


15.84 


7.89 


4 


0.42 




80424-01-01-00 


2 


2003 Aug 22.21-22.22 


0.96 


9.40 


1 


0.78 




80424-01-02-00 


2 


2003 Aug 29.58-29.61 


2.10 


14.30 


1 


0.54 




91080-02-04-00 


3 


2006 Jan 16.92-17.00 


4.92 


4.49 





0.00 


1 


91080-02-06-00 


3 


2006 Jan 17.90-17.99 


5.04 


5.21 





0.10 


3 


91080-02-14-00 


3 


2006 Jan 21.84-21.92 


4.86 


13.73 


■1 


0.50 


11 


91080-02-08-00 


3 


2006 Jan 28.63-28.71 


4.98 


7.43 


1 


0.20 


5 


91080-02-01-00 


3 


2006 Feb 04.05-04.13 


4.74 


7.94 


2 


0.86 


18 


91080-02-02-00 


3 


2006 Feb 04.52-04.55 


2.82 


5.22 





0.91 


19a 


91080-02-01-01 


3 


2006 Feb 04.57-04.60 


2.10 


5.28 





0.91 


19b 


91080-02-10-** 


3 


2006 Feb 08.38-08.46 


4.86 


25.32 


1 


0.30 


7 


91080-02-16-00 


3 


2006 Feb 11.26-11.30 


3.54 


13.50 


3 


0.60 


13 


91080-02-12-00 


3 


2006 Feb 19.18-19.22 


3.36 


20.90 


3 


0.41 


9 


91080-02-15-** 


3 


2006 Feb 20.56-20.65 


3.18 


13.51 


3 


0.55 


12 


91080-02-18-00 


3 


2006 Feb 22.00-22.04 


3.42 


13.29 


3 


0.70 


15 


91080-02-20-00 


3 


2006 Feb 23.05-23.09 


3.18 


9.82 


3 


0.81 


17 


91080-02-17-** 


3 


2006 Mar 03.26-03.40 


3.18 


8.76 


1 


0.65 


14 


91080-02-19-** 


3 


2006 Mar 04.24-04.39 


3.24 


6.90 


1 


0.76 


16 


91080-02-03-** 


3 


2006 Mar 06.20-06.29 


2.88 


4.18 





0.95 


20 


91080-02-05-** 


3 


2006 Mar 07.12-07.20 


3.30 


4.12 





0.05 


2 


91080-02-07-** 


3 


2006 Mar 08.10-08.19 


3.42 


6.05 





0.15 


4 


91080-02-09-** 


3 


2006 Mar 09.08-09.17 


3.48 


8.27 


1 


0.25 


6 


91080-02-11-** 


3 


2006 Mar 10.13-10.22 


3.42 


8.68 


2 


0.36 


8 


91080-02-13-** 


3 


2006 Mar 11.04-11.13 


4.02 


13.05 


3 


0.45 


10 



Note. — Epochs arc separated by horizontal lines in the table. 

a Obscrvation IDs ending with asterisks (**) represent two or more merged observations. 
^For 3—24 keV photons after background subtraction. 

c The number of pulse minima obtained after excluding light segments that arc consistent with a constant flux (see text). 
The orbital phase at the middle of the observation, inferred from the fit presented in § 3.1. The phases typically extend about ±0.01 in either 
direction. 

c Epoch 3 spectra number labels for Tables 4 & 5, ordered by orbital phase. Observations 91080-02-02-00 and 91080-02-01-01 were summed for 
spectral analysis. 
^This observation ended in November. 



eclipse duration implies that pulsations should be ob- 
servable at phase 0.1. Therefore, although all of the ex- 
cluded observations occurred near the eclipse, they did 
not necessarily take place during the eclipse. The light 
segments that were not consistent with a constant were 
then manually aligned to create a preliminary pulse tem- 
plate. Each light segment was cross correlated with the 
preliminary pulse template, and the cross correlation lag 
times were then used to create a more refined pulse tem- 
plate. This process was repeated until the pulse arrival 
times and the pulse template no longer changed, and was 
performed on the epoch 1 and epoch 3 data separately. 
The epoch 3 pulse template was used to find the epoch 
2 pulse minima. In total, 61 pulse arrival times were ob- 
tained. They are presented with their associated errors 
in Table 2. The final pulse templates for epochs 1 and 3 
are shown in Figure [3 It is interesting to note that the 
epoch 1 pulse profile is more peaked at <fi « 0.5-0.6 than 
the epoch 3 pulse profile. 

2.2. Error Analysis 

Statistical errors on the pulse arrival times were cal- 
culated using Monte Carlo simulations. The PCA back- 
ground errors were estimated from the unmodeled back- 



ground variations shown in Fig. 29 of Jahoda et al. 
(2006), ranging from -3-8% over 3-24 keV with 16 s 
binning. The background and errors were interpolated 
to 1 s resolution. We refer the reader to Thompson et al. 
(2006) for an explanation of the simulations. 

In addition to statistical errors, systematic uncertain- 
ties in the pulse arrival times can be caused by changes 
in the pulse profile due to varying flux levels, or varying 
levels of local absorption. To estimate these errors, we 
aligned each folded light segment with the pulse template 
using the measured pulse arrival times, and each aligned 
and folded light segment was then used as a trial pulse 
template for all of the other folded light segments. Fig- 
ure [3] shows the resulting histograms of the magnitude of 
the difference in the measured pulse minima versus the 
final templates for epochs 1 and 3. We find a system- 
atic deviation of 4.5 s for the epoch 1 data, and 5.9 s 
for the epoch 3 data. Epoch 2 consists of only two ob- 
servations, so no systematic deviation can be obtained. 
We therefore chose to apply a value of 5 s to these data. 
The systematic errors were added to the statistical errors 
in quadrature, which may slightly overestimate the total 
errors because the statistical and systematic errors are 
likely correlated. 
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TABLE 2 
Pulse Arrival Times 



Arrival Time a 


Statistical/Total 


Arrival Time a 


Statistical/Total 




Uncertainty 13 




Uncertainty 13 


00 


00 


00 


00 





2.3/5.0 


229359623 


3.1/6.7 


823 


1.6/4.8 


229910149 


2.6/6.5 


1650 


1.8/4.8 


229910977 


1.6/6.1 


5798 


1.5/4.7 


230282137 


0.5/5.9 


6627 


1.5/4.7 


230283372 


0.5/5.9 


7455 


1.7/4.8 


230284614 


0.5/5.9 


8287 


2.1/4.9 


230287911 


1.3/6.1 


11604 


0.9/4.6 


230531238 


2.8/6.5 


12439 


2.0/4.9 


230532480 


0.3/5.9 


13674 


1.0/4.6 


230533308 


1.1/6.0 


351102 


27.6/30.0 


231216181 


0.5/5.9 


353169 


11.3/12.2 


231217016 


1.0/9.2 


356875 


3.0/5.4 


231218242 


0.5/5.9 


363493 


3.0/5.4 


231334948 


2.5/6.2 


765590 


1.4/4.7 


231340752 


0.9/7.2 


771385 


1.7/4.8 


231341580 


0.7/6.8 


772214 


2.8/5.2 


231459595 


1.5/6.4 


773458 


2.7/5.2 


231460421 


1.4/6.4 


777171 


1.5/4.7 


231461665 


1.4/6.1 


/ / oUU4 


Z.o/ D.O 


zolooUoUo 


O.l/O. ( 


778837 


1.0/4.6 


231551114 


1.7/6.2 


1468308 


3.1/5.5 


231552374 


1.0/6.0 


1470385 


2.3/5.0 


232259983 


6.0/8.4 


1474508 


2.2/5.0 


232344878 


2.6/6.5 


1475756 


2.6/5.2 


232762873 


1.9/6.2 


152421254 


5.2/7.2 


232853859 


2.7/6.5 


153058868 


3.6/6.2 


232859646 


3.6/6.9 


228766823 


0.8/6.0 


232932443 


1.6/6.1 


228770532 


2.1/6.3 


232933281 


4.4/7.4 


228771378 


1.2/6.0 


232938665 


0.8/6.0 


228772626 


0.3/5.9 
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a Arrival times begin at HJD 51109.07. 

13 A systematic uncertainty of 4.5 (5.9) s was added to each epoch 1 
(3) arrival time. The epoch 2 data consist of only two minima, so no 
systematic deviation could be obtained. We chose to apply a value of 
5 s. 



3. RESULTS 

3.1. Pulse Timing Analysis 

The pulse arrival times were fitted with a seven param- 
eter model of the orbit and pulse period evolution. The 
arrival time of the nth pulse is given by 



t n = to + nP Q + — P Q P+a x sin i cos 



2ir (t n - TO) 



(1) 



where Pq is the pulse period at time ioi P is the pulse 
period derivative and is assumed to be constant, a x sini 
is the projected semi-major axis of the orbit, P or b is the 
orbital period of the system, and TO is a reference time 
corresponding to mid-eclipse for a circular orbit. The 
pulse number n is given by the nearest integer to 

tn — to tn — to 



(P) 



p a + o.5P(t n ~t Q y 



(2) 



Because the epoch 3 data has the greatest phase cov- 
erage, we first fit these pulse minima alone, yielding 
an orbital solution with a x sini = 100 ± 4 lt-s and 
P orb = 9.78 ± 0.04 days (90% confidence). Based on the 
preliminary epoch 3 fit, we attempted to pulse-connect 
the fit to the epoch 2 data from 2.4 yrs earlier, using 
the assumption of a single P between the datasets. A 



Fig. 2. — Final pulse templates for epochs 1 and 3 that were 
used to find the pulse arrival times via cross-correlation, for 3—24 
keV photons (solid line), 3-15 keV photons (dotted line), and 15-24 
keV photons (dashed line). The vertical dotted lines represent the 
phase corresponding to the pulse arrival times (Table 2). The error 
bars are lcr statistical uncertainties. 



convincing connection was made which requires no pulse 
number ambiguity over more than 180,000 cycles. The 
pulse period derivative expressed as a fractional rate is 
P/P = —7.6 x 10~ 4 yr _1 . The resulting projected semi- 
major axis remained the same, but the orbital period 



was refined to 9.74ll°:oo3 da y s - witn the additional 
constraints on the orbital period from this connection, 
we were able to connect the data through orbits to the 
epoch 1 data from 7.3 yrs before epoch 3 (272 orbital 
cycles). For this fit, however, the pulses could not be 
connected with a single P, and so the epoch 1 pulse pe- 
riod and pulse period derivative were allowed to vary in- 
dependently. The precision of the projected semi-major 
axis increased slightly with axsini = 101 ± 3 lt-s, and 
the orbital period was further refined to 9.7403 ± 0.0004 
days. Furthermore, the mass function was found to be 
f x (M) — 11.7 ± 1.2 Mq, and we define an epoch of mid- 
eclipse at HJD 53761.68 ± 0.04. The results of the fits 
are presented in Table 3, and the pulse arrival delays are 
plotted in Figure 21 The error bars in the figures rep- 
resent the 68% confidence level, and the secular change 
in the pulse periods has been removed so that the mod- 
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Pulse Minima Difference (s) 



Fig. 3. — Histogram of the magnitude of the difference in pulse 
minima versus the final template created by using all light segments 
as pulse templates for all other folded light curves for epochs 1 and 
3. 



ulation is purely due to orbital motion. By including 
an additional term in equation (1) representing the first- 
order term in a Taylor series expansion in the eccentric- 
ity (Levine et al. 2004), the eccentricity is 0.08^0 08, the 
projected semi-major axis is 96 ± 7 lt-s, and the orbital 
period is 9.7400±g;gggf days. However, the fit does not 
improve significantly, and with additional model param- 
eters the value of xt increases slightly from 1.40 to 1.42. 
The EXO 1722-363 orbit is therefore consistent with cir- 
cular, and e < 0.19 with 90% confidence. 

3.2. Orbital Modulation of Flux, Pulse Fraction, and 
Pulse Profile 

Using the orbital solution, we searched for orbital mod- 
ulation of the mean flux, the hardness ratio, and the 
pulse fraction using 20 orbital phase bins, shown in Fig- 
ure |3J This figure clearly shows the orbital modulation 
of the flux, and the general shape is consistent with the 
RXTE ASM light curve presented by Corbet et al. (2005, 
Fig. 10) and the INTEGRAL ISGRI light curve of Zu- 
rita Heras et al. (2006, Fig. 9), showing asymmetrical 
modulation of the orbital flux with a slightly higher av- 
erage flux at orbital phases when the pulsar is moving 
towards the solar system barycenter than when the pul- 
sar is moving away. The hardness ratio and pulse frac- 
tion are corrected for the Galactic ridge emission, which 
was measured from three observations occurring during 
eclipse (see below). The hardness ratio (using the PC A 
count rates in the 10-24 keV and 3-10 keV bands) is 
fairly constant except for the observations occurring near 
the eclipse, and at orbital phase 0.5. The pulse fraction 
for 0.15 < $orb < 0.85 varies between 40% and 65% in 



both energy bands, but decreases outside this range and 
is consistent with zero during eclipse. We speculate that 
the decrease in pulse fraction near eclipse ingress and 
egress may be due to smearing of the pulsations from 
Thomson scattering along the line of sight. Figure El 
shows modulation of the average pulse profile as a func- 
tion of the position of the pulsar in the orbit, created 
by folding all data in each phase bin, and accounting for 
slight changes in the pulse period and for the pulsar's 
orbital position. 

3.3. Properties of the Mass-Donating Companion 

Figure shows the inferred companion radius and 
mass, as well as the Roche lobe radius (Eggleton 1983), as 
a function of the binary inclination. The curves were de- 
rived using our orbital solution, the Corbet et al. (2005) 
eclipse half-angle measurement, and a neutron star mass 
of 1.4 Mq. As shown in the figure, the radius of the com- 
panion (i? c ) is greater than about 22 Rq, and at the limit 
of a zero mass neutron star R c > 21 Rq. The figure also 
allows us to constrain the inclination assuming the pri- 
mary underfills its Roche lobe. Since the uncertainties in 
the Roche lobe and primary radii are correlated through 
their interdependence on a^sini, simulations were used 
to find that i > 61° at the 99% confidence level. With 
this firm constraint on the inclination, the radius of the 
primary must be less than about 37 Rq and its mass 
must be less than about 22 Mq. 

The radius and mass constraints are consistent with a 
supergiant star with a spectral type between B0 and B5 
(Cox 2000). By combining photometric measurements 
by the Two Micron All-Sky Survey (2MASS) and the 
Deep Near-Infrared Survey of the Southern Sky (DE- 
NIS) in the J- and JT-bands, comparing the colors that 
are obtained with empirical (J — K)q values for B0I 
B5I stars, and applying extinction relationships (Ax/ Ay) 
from Cardelli et al. (1989), we find that the visual ex- 
tinction is Ay — 20.1 ± 0.6. The absolute visual magni- 
tude for B0-B5 supergiants is My ~ —6.5 to —6.2 (Cox 
2000), and the inferred spectral type is only consistent 
with the J and K apparent magnitudes if the distance to 
the system is between 5.3 kpc and 8.7 kpc. 

3.4. Spectral Analysis 

The epoch 3 PCA observations have excellent orbital 
phase coverage, allowing us to observe changes in the 
source spectrum as the pulsar travels through its orbit. 
For spectral analysis, all datasets occurring within a sin- 
gle 4 hr period were summed (see labels in Table 1). 
Although PCUs and 2 were on throughout the ob- 
servations, the loss of the propane layer in PCU has 
resulted in a higher background level. Since EXO 1722- 
363 is a relatively weak source, we only used data from 
PCU 2 to minimize concerns with background subtrac- 
tion. Because of the 1° field of view of the PCA and EXO 
1722-363's position in the Galactic plane (b = — 0?35), 
the spectrum consists of emission from the pulsar plus 
diffuse emission from the Galactic ridge. The Galac- 
tic ridge emission (GRE) was modeled following Valinia 
& Marshall (1998), using an absorbed Raymond-Smith 
plasma with temperature ~2-3 keV plus a power-law 
component with a photon index of ~1.8. Fortunately, 
due to the presence of eclipses in EXO 1722-363, the 
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Fig. 4. — Orbital solution to EXO 1722-363. Top and middle panels: Plots of the fits to each epoch. The epoch 1 (3) data is plotted with 
open squares (filled circles). The epoch 2 data is plotted with open triangles and have hats on the error bars. Bottom panel: All epochs 
plotted modulo the orbital period for two cycles. 



TABLE 3 
Orbital Solution to EXO 1722-363 



Parameter 


(units) 


Epoch 3 


Epochs 2 & 3 


Epochs l-3 a 


P i b 

■*■ pulse 
-fpulse 

a x sin i 

-^"orb 

fx(M) 
T0 C 

xl (dof) 


(s) 

(xitr 9 s s- 1 ) 

(lt-s) 
(d) 
(M ) 
(HJD) 


413.895 ± 0.002 
-10.66 ± 0.99 
100 ± 4 
9.78 ± 0.04 
11.2 ± 1.3 
53761.60 ± 0.09 
1.46 (28) 


419 onoo+0.0002 
41d.8933_ 0006 

-9.98 ± 0.004 

100 ± 4 
o 741 +0.005 
a - -0.003 

11.3 ± 1.4 
53761.68 ± 0.04 
1.47 (30) 


413.851 ± 0.004 413.89331°- oooe 
± 6 -9.98lg;{£S 
101 ± 3 
9.7403 ± 0.0004 

11.7 ± 1.2 
53761.68 ± 0.04 
1.40 (53) 



Note. — All errors are quoted at the 90% confidence level for a single parameter. 
a For this fit, P pu i sc and P pu i se were varied independently for epoch 1 (left) and epochs 2 &c 3 (right). 

b Pulse periods at to = HJD 51112.18662 for the fits to epochs 1, and to - HJD 53761.73126 for the fits to epochs 2 & 3. 
C T0 is defined to the epoch of mid-eclipse at $ or b — 0. 



GRE can be modeled by simultaneously fitting all ob- 
servations occurring during eclipse. As pointed out by 
Corbet et al. (2005), the epoch 1 RXTE observation 
30142-01-02-00 shows an apparent egress from eclipse 
(sec their Fig. 5). Applying our orbital solution and 
assuming a symmetric eclipse, we obtained a 1.68 ± 0.06 
day eclipse duration, consistent with 1.7 ± 0.1 days by 
Corbet et al. (2005), but inconsistent with the ^1.3 day 
measurement by Zurita Heras et al. (2006). Note that 



our uncertainty is partly due to a lack of ability to de- 
fine precisely when pulsations become apparent. Taking 
the lower limit of the Corbet et al. value at 1.6 days 
translates into 0.92 < <j> orD < 0.08. From this, we find 
that three of the epoch 3 observations took place during 
eclipse (spectra numbers 1, 2, & 20). The GRE param- 
eter values were then obtained by simultaneously fitting 
the spectra using XSPEC vll.3.2, the results of which 
are presented in Table 4. The model fit and residuals are 
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Fig. 5. — Orbital modulation of the mean flux (top panel), the 
hardness ratio (middle panel), and the pulse fraction (bottom panel) 
during epoch 3. The error bars for the top panel represent la stan- 
dard deviations (the uncertainty in the averages are less than 1%), 
and those of the bottom panel represent lcr statistical uncertain- 
ties. The hardness ratio and pulse fraction are corrected for the 
Galactic ridge emission. The pulse fraction is shown in two energy 
bands; the 3-10 keV data (solid line) have hats on the error bars, 
and the 10-24 keV data (dotted line) do not have hats on the error 
bars. 



TABLE 4 

Eclipse Spectrum for EXO 1722-363 



Parameter 


Value 


N H (xlO 22 cm" 2 ) 


7 1+ 2 ' 2 
— 2 3 


Raymond-Smith kT (keV) 


9 9+0:5 
^•^-0.3 


Normalization (xl0~ 2 ) 




Power law (PL) Index 


l-8l°o-i 


Normalization (xlO -3 ) 





Note. — All errors are quoted at the 90% confidence level for a single 
parameter. Spectra numbers 1, 2, & 20 were fit simultaneously from 
3-24 keV; x 2 1 = 118.2/136 = 0.87. The abundances of the Raymond- 
Smith plasma model were fixed at solar values. Model normalizations 
follow the definitions used in XSPEC (Arnaud 1996). 



shown in Figure 03 This model also accounts for any X- 
ray emission from the supergiant primary and scattered 
emission from the pulsar. 



With the GRE model parameters fixed at the best- 
fit values, the f 7 remaining spectra were fit with an ab- 
sorbed power law with a high energy cut-off plus a narrow 
Gaussian fixed at 6.4 keV to model iron line emission. 
Zurita Heras et al. (2006) found that phase-resolved 
spectroscopy did not show any variation of the continuum 
spectral shape with pulse phase, so pulse phase-resolved 
spectral analysis is unnecessary. Limited statistics made 
a high energy cut-off unnecessary for If of the 17 spec- 
tra; however, we believe that the cut-off is probably still 
present because it is a common characteristic of pulsar 
spectra (Mihara 1995). We therefore chose to include a 
high energy cut-off in all spectra, but allowed the cut- 
off and folding energies to vary only if an F-test showed 
that there is less than a 5% probability that the improve- 
ment to the fit due to the cut-off occurred by chance. 
With this criterion, spectra numbers 7, 9, 11-13, & 16 
had free values of E cut and £f id- In all other cases, 
the cut-off and folding energies were fixed to the best-fit 
value found by fitting all spectra simultaneously with a 
single cut-off and folding energy. The best-fit cut-off en- 
ergy was 16.2l° 2 keV and the best-fit folding energy was 
22.0l|' 7 keV. The spectral analysis results are shown in 
Table 5. As an example, we show the fit and residuals 
to spectrum number 8 (near $ orb = 0.36) in Figure El 
Also shown (see inset) are the 68% and 90% confidence 
intervals (A% 2 = 2.3 and A\ 2 = 4.6, respectively) for 
Ah versus photon index. Note that the two parameter 
confidence interval implies uncertainties in Ah (photon 
index) that are about 20% (5%) larger than those quoted 
in Table 5. After subtracting the GRE flux, we find the 
unabsorbed X-ray flux to be in the range 0.7-4.8 x 10~ 10 
ergs cm -2 s _1 from 2-24 keV. With the distance con- 
straints from above, the unabsorbed X-ray luminosity is 
L x ~ 0.2-4.4 xlO 36 ergs s" 1 over the same range. 

4. DISCUSSION OF THE SPECTRUM 

The orbital modulation of the EXO 1722-363 spec- 
trum, and especially the hydrogen column density, is in- 
triguing. Figure ^] shows Ah plotted against orbital 
phase for two cycles; a vertical shaded bar shows when 
the pulsar is eclipsed by the primary. The two observa- 
tions with the highest column densities occur just prior 
to ingress and just after egress, at a point where the pul- 
sar is viewed through the base of the stellar wind at the 
very edge of the primary. Increased absorption for these 
observations is also supported by the abrupt increase in 
the hardness ratio near eclipse (see Fig. EJ. Between 
phases 0.25-0.45, A^h is consistent at about fO 23 cm~ 2 , 
but apparently increases for the second half of the orbit. 
The error bars are somewhat large, but at phase 0.70 the 
increase in Ah is significant. It is difficult, however, to 
state whether we are observing variability in the column 
density due to differences within the orbit or differences 
between orbits. The specific orbit where each measure- 
ment was made is listed next to its corresponding data 
point in the figure. 

Phase-dependent modulation of the hydrogen column 
density has been investigated by several authors by 
studying the influence of a compact object on the stellar 
wind of the primary through observations and hydrody- 
namic simulations (e.g., Haberl ct al. 1989; Blondin et al. 
1990; Blondin et al. 1991; Blondin 1994). Stellar wind 
flowing past a compact object is affected by a number 
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Fig. 6. — Orbital modulation of the epoch 3 average pulse profile from 3—24 keV. The axis labels are omitted for clarity: The a;-axes is 
pulse phase from 0.0 to 1.0 and the y-axcs is counts s PCU~ 1 from 3 to the number in the upper left corner of each panel. The orbit 
from which the pulse profile came is listed in the upper right corner or each panel. Pulsations cannot be seen during eclipse, and can only 
be seen weakly at phases 0.1 and 0.9. The system is viewed from below, and the circle represents the orbit of the pulsar. The size of the 
projected stellar radius relative to the size of the projected orbit is based on the Corbet et al. (2005) eclipse time measurement. 



TABLE 5 

Epoch 3 Spectral Fits to EXO 1722-363 



Number 


Orb. Phase 3 


jV H 


Photon Index 


PL Norm. 


Cut-off Energy b 


Folding Encrgy b 


EW C 


xi 






(xlO 22 cm" 2 ) 




(XlO" 3 ) 


(keV) 


(keV) 


(eV) 


(dof) 


3 


0.10 






6+ 34 ' 2 
u -°-0.5 


16.2 


22.0 


< 4840 


0.70 (43) 


4 


0.15 


20tl« 


2-Ol°oI 


i 9 9 + 58.1 
J-^_6.8 


16.2 


22.0 


< 335 


0.65 (43) 


5 


0.20 




i 7 +o.4 

1 -'-0.3 


14.31^ 


16.2 


22.0 


< 170 


0.67 (43) 


6 


0.25 


sl 4 


1 «+ ' 3 
1 -°-0.2 




16.2 


22.0 


< 200 


1.14 (43) 


7 


0.30 


11 ±2 


1.2 ±0.1 


18.6+3 3 


16.6lJ; 7 


19.0t& 


174l£ 


0.68 (41) 


8 


0.36 


11 ±5 


1.8 ±0.2 


15.3l^ 8 


16.2 


22.0 


< 133 


0.71 (43) 


9 


0.41 


10 ±2 


0.9 ±0.1 


8.2+j;? 


16.511;? 


17.31, ! 4 


352l 2 4 2 


0.61 (41) 


10 


0.45 


10 ± 3 


1.4 ±0.1 


1 O Q+ 2 ' 3 
iz - y _3.9 


16.2 


22.0 


158lne 


0.93 (43) 


11 


0.50 


31 ± 14 


0.1 ±0.3 


1 0+ 2 ' 3 


12.8l°j 


14 2 +3 ' 2 


896l 2 ^ 7 


0.62 (41) 


12 


0.55 


12 ±6 


9 + ' 2 


1 8+ 1 - 6 


15.2±™ 


01 Q+8-8 
zl -°-10.2 


fi7 r + 138 
D '°-124 


1.04 (41) 


13 


0.60 


uti 


1.3 ±0.1 


10.9±|;| 


20.5l 2 ;^ 


c n +63.0 

o.u_ 4 g 


196 +82 


0.91 (41) 


14 


0.65 


99+23 
zz_ 13 


l.l+o-l 

— 0.5 


4.5+ 26 3 1 


16.2 


22.0 


< 638 


0.58 (43) 


15 


0.70 


30+? 


1.3 ±0.2 


i5.8lt;° 


16.2 


22.0 


< 202 


1.10 (43) 


16 


0.76 


16+ 14 
i0 -13 


1.5 ±0.3 


c n+2.3 
D - y -0.5 


17.6 ± 1.4 


17 0+13.2 


< 380 


0.61 (41) 


17 


0.81 


25+ 8 


1.3 ±0.2 




16.2 


22.0 


< 181 


0.82 (43) 


18 


0.86 


46+i° 2 


i 7+0.5 
L - ' -0.2 


201™ 


16.2 


22.0 


< 232 


0.85 (43) 


19 


0.91 


126tS 


i o+0.7 
1 - o -0.6 


, + 175.4 
0.0_2 2 


16.2 


22.0 


< 6490 


0.94 (43) 



Note. — XSPEC model: phabs*(raymond + power) + phabs*highecut* (power + gauss). All errors are quoted at the 90% confidence level for a 
single parameter. The Galactic ridge emission parameters were fixed at the best-fit values of the simultaneous fit to spectra numbers 1, 2, Sz 20, 
which occurred during eclipse (sec Table 4). Model normalizations follow the definitions used in XSPEC (Arnaud 1996). 

a Thc orbital phase at the middle of the observation, inferred from the fit presented in § 3.1. 

b The high energy cut-off and folding energies were allowed to vary only if an F-test showed that there is less than a 5% probability that the 
improvement to the fit due to the cut-off occurred by chance. In all other cases, the cut-off and folding energies were fixed to the best-fit value 
found by fitting all spectra simultaneously with a single cut-off and folding energy (see text). 

c Equivalent width of narrow Fe line at 6.4 keV. 
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Fig. 9. — EXO 1722-363 spectrum number 8 and residuals. Inset: 
Contour plot of Nn versus photon index, showing the 68% and 90% 
confidence intervals (A% 2 = 2.3 and A\ 2 = 4.6, respectively). 



Fig. 7. — Inclination versus companion radius and mass, plus the 
size of the Roche lobe, using the best-fit orbital parameters from 
§3.1, the Corbet et al. (2005) eclipse half-angle measurement, and 
assuming a neutron star mass of 1.4 Mq. The y-axis scale applies 
to both radius and mass although the units differ. The radius and 
mass curves have shaded error contours, and the Roche lobe curve 
has error bars. The vertical dotted line at i = 61° represents the 
lower limit to the inclination at the 99% confidence level. 
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Fig. 8. — EXO 1722-363 spectrum and residuals for the three ob- 
servations (spectra 1, 2, & 20) that occurred during eclipse. Emis- 
sion is assumed to be dominated by the GRE. The best-fit model 
was fixed for the fits to spectra numbers 3-19. 



of competing forces (gravitational, rotational, and radi- 
ation pressure), and photoionization of the stellar wind 
can create an accretion wake (Fransson & Fabian 1980). 
In addition, if the orbital separation is small enough, 
tidal distortion of the surface of the primary can lead to 
enhanced mass loss through the inner Lagrangian point 
L\ along the line of centers even if the primary does not 
fill its Roche lobe (Friend & Castor 1982). The presence 
of a tidal stream depends sensitively on the separation of 
the primary and the compact object, and only becomes 
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Fig. 10. — Best-fit hydrogen column density (Table 5) as a func- 
tion of orbital phase. The orbit from which the JVjj measurement 
was taken is labeled to the right of each point. The shaded bar 
indicates times when the pulsar is eclipsed by the primary. 



apparent when the radius of the primary is greater than 
95% of the radius of its Roche lobe (Blondin et al. 1991). 
Even if a tidal tream is not present, the simulations by 
Blondin et al. (1990) are described by the formation of 
nonsteady accretion wakes that produce dense filaments 
in the downstream wake. If the filaments are along the 
line-of-sight to the pulsar, which is likely to occur at 
orbital phases 0.5-1.0, the observed hydrogen column 
density will increase by a factor of a few to 10, qualita- 
tively consistent with the phase-dependent photoelectric 
absorption seen in EXO 1722-363. 

Models including a significant tidal stream show weak 
or altogether absent changes in the sign of the accreted 
angular momentum (Blondin et al. 1991). While the 
connection of pulses from 2003 August to 2006 March 
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using a single pulse period derivative is indicative of a 
stable accretion process over the 2.4 yr span (spinning 
up the neutron star), the orbital solution unambiguously 
shows that at least one torque reversal occurred between 
1998 November and 2003 August: The pulse period in- 
creased from 413.85 s to 414.65 s between late 1998 and 
2003 August, but has been steadily decreasing through 
early 2006. Note that the low-energy light curve of EXO 
1722-363 showed no dips, and so there is no evidence for 
a persistent accretion disk. The spin-down therefore oc- 
curred sometime between 1998 November and 2003 Au- 
gust. Such spin-down can occur via the "propeller effect" 
(Illarionov & Sunyaev 1975), or if there are significant 
azimuthal gradients in the wind density or velocity near 
the compact object (Wang 1981), however, as mentioned 
above, a torque reversal is inconsistent with the simula- 
tions predicting increases in absorption at later phases 
due to a dominant tidal stream. 

Observations of phase variability of the column density 
in Vela X-l and 4U 1700-37 qualitatively agree with that 
seen in EXO 1722-363, with significantly higher absorp- 
tion after <i> or b = 0.5 (e.g., Haberl 1989; Haberl et al. 
1989; Charles et al. 1978). A base level column density 
is seen to be reproducible over many orbits, along with a 
highly variable component causing iVjj to increase by an 
order of magnitude over the base level, usually appearing 
at later phases. This variable component is suggested by 
Blondin et al. (1991) to result from separate clumps of 
material formed by dynamical instabilities in the accre- 
tion bow shock. Interestingly, one of the narrow peaks 
in Vela X-l consistently appears at phase 0.5, and we 
see similar behavior, at least for one particular orbit, in 
EXO 1722-363. 

The picture that emerged above is one of the pulsar 
moving through the stellar wind of its B-type companion, 
affecting the distribution of matter through a number of 
competing processes, and possibly leading to a wake of 
dense material trailing the pulsar in its orbit. This ma- 
terial will also reprocess the radiation incident upon it, 
leading to fluorescent iron line emission. The interpreta- 
tion of the variability of the Fe Ka line is complicated 
because the equivalent width of the line depends on many 
factors, such as the spectrum of the continuum emission, 
and the distribution and density of the fluorescing mate- 
rial. The results of Table 5 show that the photon index 
varies over a wide range throughout the orbit (the pho- 
ton index T is consistent with values between -0.5 and 2.7 
at the 90% confidence level), as does the hydrogen col- 
umn density. The distribution of the surrounding matter 
is also expected to vary over the orbit, given the nature 
of stellar winds from massive stars and the effect of the 
pulsar on these winds. 

Using Monte Carlo simulations, Leahy & Creighton 
(1993) found that an input spectrum with T = 1 emerg- 
ing from a uniform and spherically symmetric distribu- 
tion of matter will lead to an iron line equivalent width 
of about 100 eV for a column density of 10 23 cm~ 2 . 
These authors also found that equivalent width increases 
approximately linearly below a column density of 10 25 
cm~ 2 . An input spectrum with T = 2, on the other hand, 
results in equivalent widths about one-third as large, sim- 
ply because a spectrum described by a steeper power 
law has relatively fewer photons above the iron K-edge. 
Figure shows the column density versus the equiva- 



lent width of the iron line (left panel), and the 7.1-24 
keV continuum intensity versus the iron line flux (right 
panel). Unfortunately, upper limits on the iron line in- 
tensity were obtained for more than half of the spectra 
(closed circles in the figure). Also, the data points from 
spectra 3 and 19 were not included in the figure because 
the flux for these observations is only ~10% above the 
GRE emission, and poor statistics lead to very weakly 
constrained equivalent widths. If we only consider the 
data points that are not consistent with zero (open di- 
amonds in the figure), one can see a rough correlation 
between the Fe Ka equivalent width and the Nh col- 
umn density, but the magnitude of the equivalent width 
is systematically larger than that predicted by the Leahy 
& Creighton (1993) simulations. Perhaps this reflects an 
iron abundance greater than the solar value. The depen- 
dence of the equivalent width on the iron abundance is 
linear when the optical depth is much less than one, but 
becomes less than linear when the effects of opacity be- 
come important (Matt 1997). Nevertheless, by including 
data points that are upper limits, this rough trend no 
longer exists. One factor that may contribute to the lack 
of a clear correlation is the significant variability of the 
photon index. The three largest values for the equiva- 
lent width (disregarding upper limits) are from the three 
spectra with the flattest power laws (best-fit photon in- 
dexes ofT = 0.1, 0.5, & 0.9, respectively), as we would 
expect. The right panel of Fig. ^2 shows that the cor- 
relation between the iron line flux and the continuum 
intensity above the iron K-edge is better, though still 
not excellent. If the correlation were linear and the spec- 
trum constant (minus changes in absorption), it would 
indicate that the origin of the iron line is from the same 
material as is responsible for the absorption of the con- 
tinuum. Spectral variability makes such a comparison 
impossible, but we speculate that the pulsar may be il- 
luminating matter that is not along its line of sight, and 
that may have a smaller . That the equivalent width 
increases near $ or b = 0.5 possibly indicates that a com- 
ponent of the fluorescing material is the surface of the 
primary. 

5. SUMMARY AND CONCLUSIONS 

We have found the orbital solution to the eclipsing 
HMXB EXO 1722-363 by applying pulse timing anal- 
ysis to twenty recent RXTE observations (nearly evenly 
spaced in orbital phase) and two sets of observations oc- 
curring 2.4 years and 7.3 years earlier. With the orbital 
solution, we constrained the nature of the mass-donating 
companion star, and we investigated the evolution of the 
spectrum through the pulsar's orbit. Our primary results 
are summarized as follows: 

1. The circular orbital solution is characterized by 
a x sin i = 101 ± 3 lt-s, P orb = 9.7403 ± 0.0004 days, 
and an epoch of mid-eclipse at HJD 53761.68 ± 
0.04. The eccentricity is less than 0.19 at the 90% 
confidence level. Our orbital solution also implies 
that a reversal of the torque on the neutron star oc- 
curred between 1998 November and 2003 August. 

2. Using the orbital solution and the eclipse half-angle 
measurement of 8 e — 31?8 ± 1?8 from Corbet et al. 
(2005), we constrained the companion star's mass 
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Fig. 11. — Column density versus the equivalent width of the iron line (left panel), and the 7.1—24 keV continuum intensity versus the 
iron line flux (right panel). Data points representing upper limits are indicated with closed circles, and the non- upper limit data points 
are indicated with open diamonds. The error bars in the left panel represent the 90% confidence level, while the error bars in the right 
panel represent the 68% confidence level. The data points from spectra 3 and 19, which have extremely large upper limits to the equivalent 
width, are not shown in either panel. The dotted line in the left panel shows the predicted relation between Nn and the Fe Ka equivalent 
width for power-law radiation with T = 1 emerging from a uniform sphere surrounding the source, and is the result of simulations by Leahy 
& Creighton (1993). The dashed line is the same but for Y = 2. 



and radius as a function of inclination. Assum- 
ing a neutron star mass of 1.4 M & and that the 
primary underfills its Roche lobe and, the system 
must have an inclination greater than 61° with 99% 
confidence. With the lower limit on the inclination, 
the primary's mass is between about 11 Mq and 22 
Mq, and its radius is between 21 R & amd 37 i? . 
The range of acceptable masses and radii are con- 
sistent with a supergiant star with a spectral type 
between BO and B5. 

3. There is only one 2MASS IR counterpart in the 
4" XMM-Newton error circle. Assuming this source 
is the primary, we found that the photometric mea- 
surements from 2MASS and DENIS are consistent 
with the inferred spectral type if Ay = 20.1 ± 0.6 
and the distance to the system is between 5.3 kpc 
and 8.7 kpc. In the 2-24 keV band, the unabsorbed 
X-ray luminosity is L x ~ 0.2 -4.4 xlO 36 ergs s _1 . 



4. The spectrum of EXO 1722-363 shows variable flux, 
absorbing column density, power-law photon index, 
and iron line equivalent width through the course 
of the orbit. Between orbital phases 0.25-0.45, Ah 
is consistent at about 10 23 cm~ 2 , but apparently 
increases for the second half of the orbit. 
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